Generalized Stieltjes transforms: basic aspects 

D. Karp* and E. Prilepkina^ 



Abstract. The paper surveys the basic properties of generahzed Stielt- 
jes functions including some new ones. We introduce the notion of the 
exact Stieltjes order and give a criterion of exactness, simple sufficient 
conditions and some prototypical examples. The paper includes an ap- 
pendix, where we define the left sided Riemann-Liouville and the right 
sided Kober-Erdelyi fractional integrals of measures supported on half 
axis and give inversion formulas for them. 
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1. Introduction. The generalized Stieltjes transform of a non- negative measure /x sup- 
ported on [0, oo) is defined by 



where a > and we always choose the branch of the power function which is positive on 
the positive half-axis. The measure is assumed to produce a convergent integral for each 
z G C\(— oo, 0] thus generating a function holomorphic in C\(— oo, — r], where r = inf{a; : x G 
supp(yu)}. Functions representable by the above integral plus a non- negative constant are 
known as generalized Stieltjes functions [301 132], [SSI Section 8], [3l| Chapter VIII]. 

The case a = 1 has been thoroughly studied by many authors beginning with the classi- 
cal work of Stieltjes [31] followed by Krein (see [18] and references therein), Widder [351 l3l]. 
Hirsch [TT], Berg [3] and many others. Among the most important tools facilitating such 
study are the complex inversion formula due to Stieltjes [3Tl[3l], the complex variable charac- 
terization found by Krein (see Theorem [TT] below) and the real inversion formulas by Widder 
[35],[34j. When the measure /i has compact support, the Stieltjes functions are also known 
as Markov functions studied by Chebyshev and Markov in connection with continued frac- 
tions. The deep connection of the Stieltjes and Markov functions with continued fractions 
and Fade approximation is investigated in the monographs [21 IB]- See also the survey paper 
[lOj . Connection with Bernstein functions and various other similar classes can be found in 
the carefully written recent monograph [28]. 
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For general a > much less is known. A complex inversion formula in this case has 
been found by Sumner [32] and later rediscovered by Schwarz [2S]. It was amended by 
several other complex inversion formulas by Byrne and Love in [B]. These authors also 
found several real inversion formulas in [201 121] ■ A simple real variable characterization has 
been discovered recently by Sokal [30] generalizing the corresponding result of Widder. The 
asymptotic expansion of generalized Stieltjes transforms of slowly decaying functions has 
been studied by Lopez and Ferreira in p!9], factorization as iterative Laplace transforms - 
by Yiirekli |36j, closely related classes on half plane have been investigated by Jerbashian 
|13j . An interesting connection to entire functions has been discovered in a recent work of 
Pedersen [23] . 

In this paper we collect a number of facts about generalized Stieltjes functions. Most 
of them are scattered in the literature, those we could not find are furnished with detailed 
proofs. Some of them may be new. In the last section we introduce the notion of the 
exact Stieltjes order, which is the "natural" exponent defined for each generalized Stieltjes 
function. We give a criterion of exactness and its two practical corollaries. We also added 
an appendix, where we define the left sided Riemann-Liouville and the right sided Kober- 
Erdelyi fractional integrals of measures supported on [0, oo] (the one point compactification of 
[0, oo) - see details below) and give inversion formulas for them. Our study of the generalized 
Stieltjes transforms carried out in this paper has been largely motivated by the applications 
to the theory of hypergeometric functions presented in our forthcoming work [H] . 



2. Definition and real variable properties. Define Sa, a > 0, to be the class of 
functions representable by the integral 

J [u + Z)°' 0" 
(0,00) 

where < /Uo,/ioo < oo and /i^ runs over the set of non-negative measures supported on 
[0, oo) such that 

< oo. (2j 



[0,oo) 

This condition guarantees the finiteness of the integral for all z G C\(— oo,0]. The 
classical Stieltjes cone S corresponds to the case a = 1, i.e. S = Si, see [3} formula (1)], [30l 
formula (2)] or [281 formula (2.1)]. 

If we define [0, oo] to be the one point compactification of [0, oo) formula ([1]) may be 
rewritten as 

f{z)= I {^y~^o.{du), (3) 

[0,oo] 

where 

(1 + uj" 

is a finite measure on the compact interval [0, oo]. Here Sa stands for the Dirac measure with 
mass 1 concentrated in a. This explains the notation for the non-negative constant in 
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([T]). The set of measures supported on [0, oo] and satisfying ([2]) will be denoted by A^q. The 
majority of references on the classical and generalized Stieltjes functions use formula (JT]) (or 
its particular case a = 1) to define them. See, for instance, [31 [HI El EHl ED ESI EOl EH ES]- 
However, the literature on Pade approximation frequently defines the Stieltjes functions by 
a = 1 case of the following formula 

f paidt) Poo 

(0,oo) 

Define the map on A^^ by 

[Nap]iA) := J t-''fx{dt) = J u'^fxidu) for each Borel set A C (0, oo), (5) 

l/A A 

where fl is the image measure of /i under the map t ^ 1/t and by definition [Nap]{{oo}) = 
/i({0}), [A^'o/ildO}) = /i({oo}). If (II]) and (jlj) represent the same function, the change of 
variable t = 1/u in shows that pa = N^Pa- By the same change of variable applied to 
([2]) we see that 

Pc{dt) 



, < oo, (6) 

[0,cx>) 

so that Na maps A^^ into itself. Moreover, is easily seen to be an involution on M.a'- 
NaNaP = p for each p G A^a. 

Definition (jl]) is a natural extension of the definition of Stieltjes functions used in (21 
formula(5.1)] and [TOl formula (1)]. In some situations this representation leads to simpler 
expressions for hypergeometric functions. We will work with both representations dTj) and 
(jlj). If we define the finite measure 

Pa{du) ... r 

+ PoOq + PooOoo 



'l + u) 



on the compact interval [0, oo], then pa and pa from ([3]) are related by Pa{A) = pa{l/A) for 
each Borel set Ac [0, oo]. 

We will denote by : [0, oo) — t- [0, oo) the left-continuous distribution function of the 
measure p G AIq,: -F^(a;) = p{[0,x)) normalized by F^{0) = 0. The distribution function 
defines the measure p uniquely except for a possible atom at infinity which must be specified 
separately. According to [4, section 1.8] or [271 Problem 7.9(iii)] every measure p G AIq is 
generated by such non-decreasing left- continuous function satisfying 



[0,oo) 



and a non-negative constant p^o '■= p{{oo}). Here the integral is understood as Lebesgue- 
Stieltjes integral [251 Section 20.3] or gl 2.12(vi)]. 

The Stieltjes cone S possesses a number of nice stability properties which can be found in 
[5]. The majority of these properties do not carry over to Sa, a^\. On the other hand, here 
we have some new effects related to transition from S'a to 5/3, /3 7^ a and certain stability 
properties of the class Soo '■= ^a>oSa- Below we list the basic facts about these classes. 
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Theorem 1 If f ^ Sa then g{z) := z °'f{l/z) also belongs to Sa and their representing 
measures are related by fig = N^fif. 

Proof. Indeed, using definition we get: 

Z }[1/Z)-Z j J—^-^ + Z Po + Poo- j -(—^ + Z PO + Poo 

(0,00) (0,00) 

which is precisely the representation ([T]). Since pa € Ma is arbitrary, the claim follows. □ 
The next important result is due to Sokal [30] : 

Theorem 2 (Sokal, [30]) A function f defined on (0, 00) has holomorphic extension 
f E Sa if and only if 

F„"^(x) := > (7) 

for all integers n,k > and all x > 0. Here D = d/dx. 
Remark 1. Differentiating ([1]) under the integral sign or writing 

(-l)"D^(x"+'=+"Z}"(-/'(a;)) = (8) 
we see that f G Sa implies — /' G S'^+i, so that (—!)"/(") g 6*00 for all integers n > once 

Theorem Z If a < (3 then Sa C Sp. Moreover, each f G Sa defined by can be 
written as 

(0,00) 

where pp E M.p and 

''^<*' = r(a)r(/3-a) J 

[o,y) 

Conversely, given p/3 G A^^g which represents f E Sa we can recover pa from 

r(a) f d f Pp{du) 



Y{P)V{a-P + n + l)\dy) j {y - uf—^ ^^^^ 

[o,y) 

with n = [(3 — a] and Pq({oo}) = ^^{{oo}). 

Remark 2. Formula ( [T0|) generalizes [M] Chapter VIII, Corollary 3a. 1, p. 330] which in our 
notation connects the measures pi and p2- The transformation p^ — ^ P/3 defined in (fTOj) is 
the left-sided Riemann-Liouville fractional integral. Its precise definition and inversion are 
investigated in the Appendix. 
Proof. Since 

r(/3) [ t^—^dt ^^2) 



[u + zY r(a)r(/3-a) J {z + u + ty 

(0,00) 



we have 



fia{du) r(/3) f , , , /■ t^-'^-^dt 



{u + z)» T{a)T{/3-a) J ' J {z + u + ty 

[0,00) [0,00) {0,00) 



y=u+t 



_ T{(3) f f {y-u)^—'dy 

-r(«)r(/3-a) J ^"^^"^ J {z + yy 

[0,00) (UjOo) 
(0,00) [0,?;) [0,aD) 

where iJipi^dy) is defined by ( !T0|) and we can use ToneUi's theorem [251 Chapter 20, Corollary 7] 
or [U Theorem 3.4.5] to justify the interchange of integrations. Also setting 2; = 1 we see 
by Tonelli's theorem and condition (|2]) that the measure belongs to Jvlp. According to 
Remark Al in the Appendix the measure has no atom at zero. This allows us to remove 
zero from the domain of integration in ([9]). This proves (E]) and the inclusion Sa C Sp. The 
inversion formula f lTTj) is Theorem lAll in the Appendix. □ 

Theorem 4 Each f E Sa defined by @ can also be written as 

/w= y Tf^+po. 

J {1 + XZ)P 

(0,00) 

where E M.p and 

r(/3)x"-Ma; r [ u'-^pa{du) ] ,^ ^, 

P^^''^ = Tianp-a) { J W^i)Ul^P-\- (14) 

{x,oo) 

Conversely, given pp E M.13 which represents f E So. we can recover pa from 

r(«) 



FpAy) = po + 



T{(3)T{a- (3 + n+l) 



^ dx) ^ y - x)/^-"-" ^ V ^ dy) ^ 

{0,y) {x,oo) iy,oo) 

with n = [(3 — a] and pa{{oo}) = [T{f3 — a)T{a + l)/r(/3)] lim y~°'Fp„{y). 

Remark 3. The transformation pa — ?• Pjs defined in ( IT^ is the right-sided Kober-Erdelyi 
fractional integral. Its precise definition and inversion are investigated in the Appendix. 

Proof. To demonstrate ( IT^ we employ the connection formula p*p{dy) = y~^ ppi^dy), 
where is the image of under y — )• y^^. We have {A = T {(3) / [T {a)T {f3 — a)] for brevity) 

y^p;idy) = ppidy) = Ady J .^J^^^L_ + p,Ay^'-~^dy 

(o,y) 

(l/?/,oo) (l/j/,oo) 



(15) 
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Dividing by and making substitution y = 1/x, dy = —dx/x^ we arrive at (JH]). The 
inversion formula f llSp is Theorem IA2I in the Appendix. □ 

The following result is also due to Sokal [SDl formulas (10a) and (10b)]. 

Theorem 5 na>o ~ {non-negative constants}. 

This result suggests the following definition: 



Recall that a function / : (0, oo) — )■ M is said to be completely monotonic if / has 
derivatives of all orders and satisfies (— l)"/*^"')(x) > 0, for all x > and n = 0, 1, . . . We 
denote the set of completely monotonic functions by CA4. The following result was pointed 
out to us by Christian Berg and is also hinted at in [30] . 

Theorem 6 Soo = CM., where the closure is taken with respect to pointwise convergence 
on (0, oo). 

Proof. The inclusion S'oo C CA^ follows from the fact that each / G Ua>o completely 
monotonic combined with closedness of CAi under pointwise convergence [281 Corollary 1.6]. 
To prove the reverse inclusion recall that according to Bernstein's theorem (see [3U Chapter 
IV, Theorem 12b] or [28l Theorem 1.4]) each completely monotonic function is the Laplace 
transform of a nonnegative measure: 



Since f{x) is non-increasing the sequence /(1/n) is non- decreasing. Hence, two cases are 
possible: 1) /(1/n) is bounded by a constant or 2) lim„_j,oo f{^/n) = +oo. Define a„ := y/n 
in the first case and a„ := /(1/n) in the second. Define the following sequence of functions: 



5*0 := {non- negative constants}. 




[0,oo) 




[0,oo) 



Obviously, 




[0,cx.) 



SO that the integral defining /„(x) exists. Moreover, 




[0,oo) 



where an{du) = af^e~^a{du). Clearly, /„ G Ua>o'^"- Next, we have 




— e 



—xt 



+ e 



—xt 



— e 
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f^{x)-f{x)= I + '^'^ -e"Ae'^a{dt) + f{x + l/n)-f{x). (16) 



[0,00) 

It is verified by straightforward calculus that the maximum of 



m ■- 

is attained at the point t* satisfying 



rf \ """"^ 



Hence, 

max|V'(t)| = ^ < — 



It follows from (fT6|) that 

|/n(x) - < + \f{x + 1/n) - 



Due to the definition of a„ and continuity of /(x) we conclude that lim„^oo fn{x) = f{x) for 
each a; > 0. □ 

Theorem 7 If f ^ Sa and g E Sp then fg G Sa+i3- 

Proof. See [121 Chapter VII, paragraph 7.4]. 

Remark 4. Theorems [3] and [7] show that the union 6*00 is a cone with multiplication: if 
f,g^ Soo then af + bg E Soo, for all a,b > and fg G 5*00. 

Theorem 8 Each 

f{z)= j (^7^ + ^^e^« 

[0,00) 

can be represented in the form 

f{z) = u)du; z) + /ioo, (17) 

T[a) 

where C denotes the Laplace transform. 
Proof. Write 

1 1 



{u + zY r(a) 

[0,00) 



and apply Tonelli's theorem to show that the iterated integral in (|T71) exists and is equal to 
Remark 5. Formula ( !T7|1 has been found in ^36j for absolutely continuous measures. 
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Introduce the standard notation for the right-sided Riemann-Liouville fractional integral 
[T5| section 2.2], [26i §5] and the right-sided Caputo fractional derivative [151 section 2.4], 
section 2.4.1]: 

oo 

f{t)dt 



r(A) J it-x) 



l-A • 



-D-f = (n-V = ^^ /-^^^^^ n=[A] + l 

a; 

Theorem 9 For a fixed non-negative measure ^ G Aia CLnd (3 > a > denote 



[0,oo) 



and 



Then for all x > 



[u + z 

[0,oo) 



+ Hoc- (19) 



r(/3) /• (^(t) - ^oo)dt , r(/3) 
= r(c.)r(/3-a) y (t-x)^w + = VJ^(^) - /^oo](x) + /.oo (20) 



and 



f«(x) - / /^^(^)^^ + . - mcD- If m,) + a (21) 



r(/3)r(n - /3 + a) 7 (t - x)i+/5-"-" '^"^ r(/3) 

where n = [f3 — a] + 1. 

Proof. To prove (1201) substitute ( lT9ll for //3 into (l20l) and exchange the order of integration 
which is legitimate by Tonelli's theorem again. Conditions P > a > guarantee the existence 
of the inner integral. Formula ( 12T|) is one of several inversion formulas for the Riemann- 
Liouville integral in fl20|) which is applicable since fa is infinitely differentiable (see, for 
instance, [151 section 2.4], [211 section 2.4.1]). To demonstrate its validity differentiate under 
integral sign (n > 1^ 



fi-\t) = (-l)"(a)„ 



fi{du) 



(M + t)"+'^ 
[0,oo) 

where (a)„ = a{a + 1) ■ ■ ■ {a + n — 1) = T{a + n)/T{a), substitute into f l2T|) and exchange 
the order of integrations. □ 

Remark 6. The Riemann-Liouville fractional derivative cannot be used in (1211) since, in 
general, the resulting integral would diverge. 

Theorem 10 The class Sa, a > 1 is closed under pointwise limits: if {fn}'^=i C Sa and 
if the limit lim fn{x) = f{x) exists for all x > then / G S'q. 

n— >oo 

Proof. A proof for a = 1 can be found in [ill Proposition 1] or by a different argument in 
[281 Theorem 2.2(iii)]. The latter carries over mutatis mutandis to all a > 0. □ 
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3. Complex variable properties. Clearly, / is holomorphic in C\(— oo,— r], where 
r = inf{x : x G supp(/iQ,)} and f{^) = f\z). In particular, if ^ supp(/iQ,) (or, equivalently, 
supp(pq) is bounded) then the function / can be represented by the power series 

/(.) = Y.^-ir^Pu{a)z\ 

k=0 

convergent in the disk \z\ < 1/R. Here R = sup{x : x G supp(pa)}, {a)k = r(a + k)/T{a), 
and 

Pkia)= j t'^dp^it) <oo, A; = 0,1,... (22) 

lo,R] 

are the moments of the measure pa which are finite due to ([6]). 

For functions belonging to Si Krein [T8| Appendix, Theorem 4] found the following 
celebrated characterization. 

Theorem 11 A function f holomorphic in the cut plane C \(— oo,0] belongs to Si iff 
f{x) > for X > and '^f{z) < for '^z > 0. 

Because of the special role played by Si (including a number of stability properties [3], 
connection to continued fractions [8] and Fade approximation [2]) it is interesting to relate 
the functions from Sa to ^i. One way of doing this is provided by Theorem |9] (it works if 
/i G Ml), another approach is presented in the following two theorems. 

Theorem 12 Suppose f e Sa, < a <1. Then /^Z" G Si. 

Proof. For '<sz > and t > we have (arg(z) denotes the principal value of the argument 
of z): 

-na < arg{z + t^ < ^ - ixa < aig{f{z)) < ^ $5(/(2)^/") < 0. 

Since f^^"{z) is holomorphic in C\ (— oo, 0] and non-negative for z > we get the conclusion 
by Krein 's theorem [TTl □ 

Theorem 13 Suppose f e S^, a>l. Then g{z) := /(-z^/") G ^i. 

Proof. Indeed g{z) is holomorphic in C \ (— oo, 0] and non- negative for x > 0. Next for 
'^z > and t > we have (arg(2;) denotes the principal value of the argument of z): 

< arg(;z^/") < vr/a < arg(2^/" + t) < n/a, 5>(;z^/" + 1)"" < 

Integrating the last expression with respect to non-negative measure preserves the lower half 
plane, so that the proof is completed by Krein's theorem [TTl D 

Remark 7. For a > 1 the mapping Sa Si defined by f(z) g(z) := f{z^^°') is clearly 
not surjective as can be seen immediately by taking g{z) = 1/(1-1-2;) G Si, f{z) = gi^z"^) ^ Sa 
(since it is not holomorphic in the upper half-plane). The conditions (a) f{z) is holomorphic 
in the cut plane C\ (— oo, 0], (b) f{x) > for x > and (c) ^f{z) < for < arg(z) < n/a 
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are necessary for / to belong to Sa- Unfortunately, these conditions are not sufficient as the 
following example shows: 

l__ 

{z + iy 2(z + 2)2- 

Indeed, a straightforward computation yields {z = x + iy)\ 

?/(30 + 87x + 96x2 ^ 5Q^3 ^ i2x^ + + 12y^ + 22xy^ + Ux'^y^ + 2x^y^ + xy^) 



Hence, for {z E C : ^z > 0,'^z > 0} we get '^f{z) < and f{x) > for x > 0. However, 
f ^ S2 since the representing measure is signed. In fact, f E S3 since 



where du(t) = l^^'^ '^^^'^l' 
wnereaz/^rj \ dt, t G (2, 00). 



We thank Alex Gomilko for this example. It provides a partial answer to the question asked 
by Sokal at the end of [30] . 



Theorem 14 Each f E Sa satisfies 

\m\<A 



z - 1 



Q(z) 



(23) 



where A = L . fia{du)/{u + 1)" < 00, and ^{z) < 0. 



'[0,00) 
Proof. We have 



1/(^)1 



[0,00) 



j2a{du) 

(u + zY 



fla{du) fu + 1 



[0,00) 



u + 1 



< A max 

u>0 



u + z 



+ /ioo = A[7/;(^)]"/2 + /i, 



with A given above and 



z I = max 

ii>0 



u + 1 


2 


z-l 


u + z 




%z) 



The last equality is true for '^{z) < by standard calculus. □ 

Remark 8. The estimate fl23|) may seem very weak because the right hand side does not 
tend to zero as |z| — 00 while f{z) — )■ 0. However, the decrease of / may be arbitrarily slow 
so that it is difficult to expect a better estimate valid for the whole class S^- 
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4. Exact Stieltjes order. We will say that / is of exact Stieltjes order a* if / G 6*00 and 

«*[/] = inf{a: / G 5J. (24) 

Theorem 15 /// is of exact Stieltjes order a* then f G S^*- In other words the infimum 
in ( 124|) is always attained. 

Proof. Since / G 6*00 it is infinitely differentiable so that F^f,{x) in ([7]) is well defined and 
continuous in a for each fixed x > for all a > a*. Passing to the limit a — )■ a* we verify 
that ([7]) is true for a = a*. Hence, by Sokal's theorem |2]/ G Sa*- □ 

Theorem 16 Suppose f & Bp. Then /3 > «*[/] iff the function 

(o,y) 

is non- decreasing on (0, 00) for some e G (0,min{/3, 1}). 

Proof. Assume $(?/) is non- decreasing. We need to show that /3 is not exact. According 
to Definition 1 in the Appendix and Theorem lAll 



^{jj)dy 



where a = {3 — e and /loo = /^/^({oo}). Hence, 



[0,00) 

This implies that 



X , 1 < C)0. 



<oo. (26) 



;i+i/)" 



[0,00) 

Indeed, for each t > integration by parts, 

$(7/) * 1 f d^{y) f ^{y)dy 



a{l + yy 



^ a J {l+yY J (1 + 1/)"+^ 



shows that limt_^oo '^'(^)(1 + t) " exists (finite or infinite). This limit must be zero since 
otherwise $(|/)(1 + y)"" > C > for all ?/ > M and 

> C / = 00. 



(1 + y)"+i J 1 

[0,00) [M,oo) 

Hence, limt_>oo $(t)(l + 1)~° = and the above integration by parts proves fl26|) . It follows 
that the measure fia whose distribution function is equal to A^[y) and whose atom at infinity 
is equal to fi^o belongs to A^^. Here A = T{a)/[T{(3)T{a — (3 + 1)]. Consider the function 

'^'^= J (^7T^ + ^- 

[0,00) 
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By Theorem [3] we have g E and 

[0,00) 

where jliji^du) is given by ( ITOj) . But then /i^ and are related by (ITT]) which coincides with 
(!25|) times A (note that n = in ( ITTj) because (3 — a = e < 1). This proves that jlp = yU^ so 
that f{z) = g{z) G S'^, a < (3. 

Conversely, if /3 is not exact choose e E (0, /3 — a*), e < 1. We have / G S'/3_e. According 
to Theorem |3] the function A^{y), where $ is defined in (l25ll equals the distribution function 
of the measure fips and so is non-decreasing. □ 

Corollary 1 Suppose f G 5*^, e G (0,min{/3, 1}) and the following limit exists: 

where $ defined by fl25l) . If A < 1 then (3 is the exact Stieltjes order of f . 

Proof. Clearly, the condition A < 1 implies that cannot be non-decreasing so that by 
Theorem [T6] must be exact. □ 

Corollary 2 Suppose / G 5/3 and the support of the measure iip is compact. Then (3 is 
the exact Stieltjes order of f . 

Proof. Indeed, for all y > B, B := sup{a; : x G supp{dfi)}, the function $(?/) is strictly 
decreasing for each e G (0,min{/3, 1}), so that by Theorem [T6l 3 must be exact. □ 

Consider three prototypical examples. 

Example 1. Find the exact Stieltjes order of {a > 1) 



dt 



Method I: by corollary |2] a* = a since supp((i/i) is compact. 

Method II: by theorem [T6] compute $(?/). Let I{A) be the indicator function of a set A. 
We have (0 < £ < 1) 



' li[0,l])du ^ fy-V(l -e), 0<y<l, 



It is straightforward to check that this function is decreasing for y > 1 so that again a* = a. 
Example 2. Find the exact Stieltjes order of {a > 1) 



dt 



(z + ty 
12 



By theorem [16] compute 

' I{[l,oo))du (0, 0<t<l 



{t-uy i(t_ 1)1-7(1 t>i 



This function is non- decreasing on [0, cxd), so that a* < a. To find the exact order we 
compute 

oo 



(a- l)(l + 2)°-i J {z + ty-^' 



where the measure dv{t) is concentrated at one point t = 1 with z/({l}) = (a — 1)^^ so that 
by Corollary [2] a* = a - I. 

Example 3. According to Euler's integral representation [H Theorem 2.2.1] the Gauss 
hypergeometric functions can be written as 

Tic) u^'^il - uY'^^^ 
f{z) :=2Fi{a,b]C]-z) = —— -/ ' ^ ' du, c>b>0. 

Assume that < a < 6. By the above formula f E Sa- Change of variable u = 1/t yields 

^/ , , r(c) f°° fi(t)dt 

^^^^"'''^'-^) = rxw^/ (JT^' 

where 

n{t) = t^-%t - iy-''~\ t > 1. 

We aim to show that a is the exact Stieltjes order of /. Computation gives 

y 

*^^) ■= I = - -a,c-b;c-b-e + l;-{y-l)), 

1 

where A = r(c — b)r{l — e)/[r{c — b — e + 1)] for any given e G (0, min{l, c — b}). Using the 
asjTiiptotic formula [U formula (2.3.12)] for 2-^1 as y — )■ 00, we get 

^y)=Aiy-l)-' P + o(l)), 

where D is a constant. This implies that limj^_j.oo[$(2?/)/$(y)] = 2^^ < 1 so that by Corol- 
lary [T] the exact order is equal to a. 

In our forthcoming paper [1A\ we will investigate the exact Stieltjes order of the general- 
ized hypergeometric function q+iFg. 

Remark 9. According to Theorem [10] the class Sa is closed under pointwise limits. It is 
then reasonable to ask whether the exact Stieltjes order is also preserved by such limits. The 
following example shows that the answer is negative in general. Consider the sequence of 
functions {a > 1) 

m 

= / = ^a-l){z + l)-i ~ {a-l)iz + m)"-i ' "^ = 2'3,--- 

1 
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According to Corollary |2] each fm has exact order a while pointwise 

1 

lim fmiz) = -, TT ^ — r for each z G C\(— oo.Ol. 

According, to Example 2 above, the limit function is of exact order a — 1. 
Remark 10. Theorem [9] shows that if a is the exact Stieltjes order of / then its fractional 
derivative of order 7 will have the exact order a + 7 while its fractional integral of order 7 
will have the exact order a — 7 provided that a > 7 and [i G A^a_-y. 
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APPENDIX 

Inversion formulas for some fractional integrals of 

measures on half-axis 



In this appendix we prove some facts about fractional integrals and derivatives of Borel 
measures supported on IR+. Fractional calculus is certainly a classical subject with a number 
of great monographs available today, including [9l [151 HSl 1221 [Ml [26] . Fractional integrals and 
derivatives have been studied in (weighted) spaces of integrable functions, in Holder classes, 
in spaces of generalized functions, in the complex plane, for functions of several variables 
and in various other contexts. We could not find, however, a good reference for fractional 
integrodifferentiation of measures. One reason could be is that the definition of the fractional 
integral of a measure is fairly straightforward. The inversion problem, nonetheless, i.e. the 
definition of a fractional derivative, might not be so trivial. Here we prove two inversion 
theorems required in the study of generalized Stieltjes transforms. They might be useful in 
some other contexts as well, for instance in connection with completely monotonic functions 
of positive order, see [T7] . 

Let us remind the reader that A^q is the positive cone comprising non-negative Borel 
measures supported on [0, oo] and satisfying (|2|). Clearly, A^aCA^/3ifa</3. There is a 
natural involution Na defined on A^q, by (|5|). For a measure /i G A^q, we denote by its 
left-continuous distribution function normalized by F^{0) = 0. The distribution function 
defines the measure uniquely except for a possible atom at infinity which must be specified 
separately. 

Definition Al. Let fi G Ma- The measure u := J^yU is called the left-sided Riemann- 
Liouville fractional integral of fi of order 77 > if 

u(B) := -i- ! dy f , + ^[B n {00}) for each Borel set B C [0, oo\. (27) 

^wi J J [y-uy 

B\{oo} %y) 

Remark Al. Formula fl27p is certainly a straightforward generalization of the left-sided 
Riemann-Liouville fractional integral as given in jl5i (2. 2.1), (2. 2. 2)] and [2B], Chapter 2(5.1), (5. 3)]. 
One can check that z^({0}) = regardless of /i({0}) by computing the limit limj^_j.o Fi^iy) = 0. 

Definition A2. Let fi G Aia- The measure r := K~^fi is called the right-sided Kober- 
Erdelyi fractional integral of /i of order 77 > if 

riB)■.= ^^iBn{0}) + -^ [y^-'dyl / — -^^^^^ + k (28) 



^iv) J ] J u'l+'^-^iu-yy' 

B\{0} [iy,oc) 

for each Borel set Be [0, 00]. 

Remark A2. Formula ( |28l) is certainly a straightforward generalization of the right-sided 
Kober-Erdelyi fractional integral as given in [151 (2.6.8)] and |26l Chapter 4, (18.6)] with a 
slight change of notation. Condition (|2|) ensures that the right hand side of ([28]) exists. By 
definition the measure r has no atom at infinity. 
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Theorem Al ///U G Aia then v = I^^ G M.a+q CLiT'd given v we can recover fi from 

^"^^^^Til + n-v) (rf^) / (y-n)t"' ^^^^ 

where n = [q] and yu({oo}) = z/({oo}). 
Proof. To show that u G M.a+r] compute 

v{dt) 1 f dt f fi{du) 



[0,00) [0,00) [0,t) 

[0,00) {«,cxd) [0,00) 

V ' 

={l+M)-"B(j?,a) 

according to ([2]). Here 

r(c.)r(r/) 

B[r],a) = — — ■ — - 
T{a + rj) 

is Euler's beta function. The interchange of the order of integrations is legitimate by Tonelh's 
theorem [271 Theorem 13.8], [251 Chapter 20, Corollary 7] or [H Theorem 3.4.5]. The proof of 
(!29|) is a paraphrase of the standard proof of the inversion formula for the Riemann-Liouville 
fractional integral [261 Theorem 2.4] except that we recover the distribution function of 
the measure so that we differentiate one time less the standard and we employ Tonelli's 
theorem to justify the interchange of integrations. Substitute fl271) into the integral on the 
right hand side of fl29|) (recall that n = 

v{du) 



[0,S/) 



1 [ du f ii(dt) 1 [ , , s. [ du 



^iv) J iy-uyi-'' J {u-ty-'i T{r]) J J {u-ty-'^iiy - uyi-"^ 

[0,j/) [0,m) [Q,y) {t,yh) 



r{l — rj + n) 



j {y-t^fiidt) =T{l-ri + n) jdt^ jdt^--- j F^{tn)dtn. 



I0,y) [0,y) [0,ti) [0,t„_i) 

Since F^(tn) is non- decreasing it is locally Lebesgue integrable which implies that the function 
on the right belongs to y4C"'[0,i?] (the function and n — 1 its derivative are absolutely 
continuous) for any R> 0. Hence, we can recover by n-fold differentiation. □ 

Theorem A2 If ^ E M.a then t := K~^{^) G M.a+r] o-nd given t we can recover fi from 

1 



F,{y) = rim 



r{n + 1 — rj) 



J \ dx J J s°+"(s - x)*?-" \ dy J j s'^+'^{s - yy-"^ 

(0,j/) (x,oo) {y,oo) 



(30) 
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where n = [rj] and fi^o = Ci^iv) 1™ V "Prill)- 

J/— >oo 

Proof. To show that r G M.a^r^ compute 

T{dt) 1 [ f^-^dt f i2{du) 



(0,00) [0,00) {t,oo) 

(0,00) (0,u) (0,00) 

V ' 

=u^+a-^{l+u)-" B{a,ri) 

The interchange of integrations is again justified by Tonelh's theorem. To prove formula 
(15U]) assume for the moment that /ioo = and substitute into (|5UI) : 

."-1^^/ ^^^Y^V^nf ^jds) f ody f T{ds) 



' ^ dx J J s"+"(s - x)*?"" \ dy J J - y)''- 



{0,y) (x,oo) (sjOo) 



T{r]) V ay/ J s°'+'^{s-yY-'' J - s)!"'' 

(y,oo) (s,oo) 
« [ ^a^l^J f dV u^V-n [ ds f fxidu) 



T{ri) J [V^V J s"+i(s - l/if:)'?-" J - s)!"'' J, 

(0,S/) (l/i,oo) (s,oo) 



t=l/x 



y"" ( / dy _ f ds f fi{du) 



t=i/j/ 



(l/t,oo) (SlOo) 

where we have used the formula 

Further, exchange of the order of integrations justified by Tonelli's theorem yields: 
ds f fi{du) 



\t=l/y 



(l/t,oo) (•5,00) 



(l/t,oo) (l/t,oo) 



=t';-"«'J-"-i(ut-l)"_B(»7,l-»;+n) 



We will show now that 

dy f fi{du 



dtj J ^ ^ 7 

{l/t,oo) (l/t,oo) 

17 



Denote by /i* the image of the measure /z under the mapping X{u) = 1/u, so that for each 
Borel set A C (0, oo) we have := /i(A-^(A)). Then 



n 



(l/t,oo) (O.t) (0,t) 

"^'(^) /^^^ jdt2--- J s''fi*{ds) = n\ J s''fi*{ds) = n\ J u-"^i{du) 

(O.t) (O.ii) (0,t„) (0,t) (l/t.oo) 

Finally, by integration by parts for Lebesgue-Stieltjes integral (see [3 Theorem 6.2.2]) we 
obtain: 

J \ dx J J s°+"(s - x)*?-" \ dy) J s"+"(s - 

(0,5/) {a;,oo) (y,oo) 



(0,2/) {a;,oo) (»/,oo) 



r(l-77 + ri) x" u-''n{du) 



u "^{du] 



+ y jj,{dx)—y"J u '^jj,{du) 

(x,oo) (0,y) (y,oo) 

= r(l-r7 + n) y" ^(rfx) = r(l-r/ + n)(F^(|/)-F^(0+)) = r(l-r7 + n)(F^(l/)-/i({0})). 

(o,y) 

Since /i({0}) = r({0}) by ( l28l) this proves formula (l30l) . While integrating by parts we have 
also used the following limit: 

limx" y"M->(rfM) = lirnx" y"M-"rfF^(M) = limx"(^M~"F^(M) +a j u-''-^F^{u)du^ 

{x,oo) (3^,1) (^,1) 

= limfax° / u--~'F,{u)du - F,{x)] = lim ^ " - F,(0+) = 0. 



Here the first equality is due to ([2]), the second is integration by parts and the preultimate 
is L'Hopital's rule applied if f^^ ^^u~^~^F^{u)du is unbounded. 

If /ioo 7^ we need to add the following term to the third line of flMl) : 



-B{T],l—ri+n)y 



r(r7) 7 \ dx J J s"+"(s - x)^-" r(r7) \ dy J J - y)''- 

(0,1/) {x,oo) iy,oo) 

" v ' 

=B{ri,l-ri+n)x-'i 

(0,y) (0,y) 



18 



where 

d \ „ \ I a \ 



dx J \ \dt 



-x^— 1 = U — ] r)- =n\. 

t=l/x 



This shows that (!3T|) is still valid for /ioo 7^ 0. 

Finally, to recover the atom at infinity /ioo we compute 



hm x-^F^x) = hm [/i({0})a;-"] + -i- lim x"" /y^-^y / — ^ 



n{du) 



(0,x) (j/,00) 



~'~ r(?7) a;-^oo ^ a;r(r7) ~'~ r(r7) i^-^oo ^ J u'^+'^~^ {u — y)^~^ 

In order to show that the last limit is zero we interchange the order of integrations (justified 
again by Tonelli's theorem): 



(0,x) {y,oo) (0,x) (0,m) 



+ J u'-^-yidu) J {u-yr-V-'dy =F^{x)-i,{{Q}) 



[1,00) (0,x) 



={l/a)x"u'^-^2Fi{a,l-r];l+a;x/u) 
1 /■ _ 

a J 

[XyOO) 

where 2F1 is the Gauss hypergeometric function [H Chapter 2]. Hence we need to prove 
that: 

lim x-'^FJx) = 

3;— >oo 

and 

lim / u~'^2Fi{a,l — rj;! + a;x/u)fi{du) = 0. 
x—^00 J 

[a;, 00) 

Both equalities follow from ([2]): the first was proved by Widder [311 Corollary 3a. 3], the 
second follows from the fact that 2-^1 (a, b; c; x) is bounded on [0, 1] if c > a + 6 by the Gauss 
formula fl] Theorem 2.2.2]. □ 

Remark A3. The operator {—x'^Dx)^ encountered in fl30|) can be expanded as follows 

n 

{-x^DxYf = ^a(n,m)x"+'"/(™)(x), 

m=l 

where the numbers 

a{n,m) = {-iy—[ 

ml \m — 1 



19 



are known as Lah numbers [3? , A008297] satisfying 

a{n + l,m) = {n + m)a{n, m) + a(n, m — 1). 

Applying —x^D^ to the above expansion we see the same recurrence which given the same 
initial values furnishes a proof of the above expansion. 

Remark A4. Another way to obtain a representation for yU via r for r/ > 1 is the following. 
Denote jl = K~i{^), i.e. according to fl28|) 



= /i({0})(5o + ^dx J u °'fi{du) = fi{{0})6o + (j)ji{x)dx, 

(x,oo) 

where the last equality is the definition of (pfxix). It is easy to verify that this formula is 
inverted as follows 



- 0) - F^(0+) = a J jl{dx)-yMy)- 



{0,y) 

If 7] > 1 we have according to fl28l) for p, G M.a+i and the semigroup property of Kober- 
Erdeliy operator [15l (2.6.24)] 

T{dx) := K-^in) = 

(x,oo) 

Here we have the standard Riemann-Liouvile fractional integral of the function m^^"^''</)^(u). 
We cannot, however, use the Riemann-Liouvile fractional derivative to invert the above 
formula, since, in general the integral in its definition will diverge. Instead, we can employ 
Caputo's fractional derivative [151 section 2.4] to invert (!32|) : 

U '(J)r.(u)= D' \X (Pr(x]\{u] = — / -^^ z , 

(m,oo) 

where n = [r]] and 

^.^ . (-1)"^°+"-!^^ f [x-"0,(s)]('^)da; 

ix{du) = ixi{0})5o ' I I ^ \ Ji 

1 i^/t — 7/ -r -l; 

(ii,oo) 

so that 



T{n-r] + l) J (x-m)'?-" ' 



y = {0} 5o + — T u du / ^zz^dx 

i[n — 7] + 1) J J [x — uyi " 

{0,y) {n,oo) 

(_l^nya+r, |. ^^-a^^ (^^^ j (n) 



r(?2 — ?7 + 1) J {x — uY' 

This formula is, however, less general then (130|) since it requires that 0r(x) G y4C"(0,oo) 
which is not guaranteed by /x G A^o. For < ?7 < 1 both formulas take the same form. 
The relation between and is revealed in the following theorem. 
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Theorem A3 Suppose ji G Ma- Then N^+r^I^iJ, = K^.^Nafi. 
Proof. Indeed, if / is given by ([1]) and v = ^ then by Theorem [3] 

jj,{du) r(a + 77) f u{du) r(a + ?7) f Ti{dt) 



. _ f fi{duj _r{a + r]) r u{duj _r(a + ?7) C 



[0,00) (0,00) (0,00) 

where ti = N^+rji^ = Na+rjI^fJ'- On the other hand, if p = Nafi then according to the 
comment after formula ([5]) and by Theorem H] we have 

_ r p{dt) HQ r(a + 77) f T2{dt) 



{i + tzY 2" r(a) J {i + tzY+'^' 

(0,00) (0,00) 

where T2 = K~.^p = K~,^Nafi. Comparing these two formulas we conclude that ri = T2 due 
to uniqueness of the representing measure. □ 
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